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Abstract 



We consider a D-dimcnsional cosmological model describing an evolution of Ricci-flat factor spaces, Mi, . . . , 
M n (n > 3), in the presence of an m-component perfect fluid source (n — 1 > m > 2). We find characteristic 
vectors, related to the matter constants in the barotropic equations of state for fluid components of all factor 
spaces. We show that, in the case where we can interpret these vectors as the root vectors of a Lie algebra 
of Cartan type A m — slim + 1, C), the model reduces to the classical open m-body Toda chain. Using an 
elegant technique by Anderson Q for solving this system, we integrate the Einstein equations for the model 
' and present the metric in a Kasner-like form. 

00 ! PACS numbers: 04.20.J, 04.60.+n, 03.65.Ge 
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a" 1 Introduction 

Recently, investigations on multidimensional gravitation and cosmology have found renewed interest. It 
is well known now, that multi-scalar-tensor models derived from a higher dimensional multidimensional 
Einstein action are similar to the (bosonic sector of) effective low-energy models from (super) string theory. 
Beyond that fact, it was recently shown that the multidimensional ansatz provides also a natural clue to 
■ membrane theory (as the natural generalization of string theory). 

Here however, we restrict to multidimensional (spacially homogeneous) cosmology 29 . The D- 



dimcnsional cosmological model describing the evolution of n (n > 3) Ricci-flat spaces Mi, . . . ,M n in the 
presence of an m-component (n — 1 > m > 2) perfect fluid source is considered. The barotropic equations 
of state for the mass-energy densities and pressures of the components are given for each space. When the 
vectors related to the constants in the barotropic equations of state can be interpreted as root vectors of the 
Lie algebra A m = sl(m + 1, C), the model reduces to the classical open-chain m-body Toda system. Using 
the new elegant form of its exact solution proposed in , we integrate the Einstein equations for the model 
and present the metric in the Kasner-like form. 

2 The model and the equations of motion 

The D-dimensional space-time manifold M may be the product of the time axis R and n manifolds 
M l7 ...,M„, i.e. 

M = RxM t x...xM„. (2.1) 
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The product of some of the manifolds Mi, . . . , Mt, 1 < k < 3, gives the external 3-dimensional space and 
the remaining part M^+i, ■ ■ ■ , M n stands for so-called internal spaces. Further, for sake of generality, we 
admit that dimensions Ni — dimMj for i = 1, . . . , n are arbitrary. 
The manifold M is equipped with the metric 

n 

ds 2 = -c^^dt 2 + cxp[2x'(t)}dsl (2.2) 

i=l 

where j(t) is an arbitrary function determining the time t and ds 2 is the metric on the manifold Mj. We 
assume that the manifolds Mi, . . . , M n are Ricci-flat, i.e. the components of the Ricci tensor for the metrics 
ds 2 are zero. Under this assumption the Ricci tensor for the metric (2.2) has following non-zero components 



- ^ N^f + 7o - 77oJ (2.3) 

iC/ = e-^[x i + x i (i -j)]6™; ( 2 - 4 ) 

with the definition 

n 

7o = 5^x\ (2.5) 

i=l 

Indices m t and n, in (2.3), (2.4) for i = 1, . . . , n run from (D - £" =i TV,) to (L» - £" =i ^ + ^i) (- = 
1 + = dimM). 

The source of gravity shall be a multicomponent perfect fluid. The energy-momentum tensor of such 
source comoving coordinates reads 

m 

2# = (2.6) 

(7ff W ) = diag -^>(t), p< s) (t), . . . . . . , p^(t), . . . , (2.7) 

\ Nitimes N n times J 

Furthermore we suppose that for any s-th component of the perfect fluid the barotropic equation of state is 
given by 

p<f\t)= (l-ft«)p«(t), s = l,...,m, (2.8) 

where h\ s ' = const. 

The equation of motion X/mT^ 1 ^ = for the s-th component of the perfect fluid described by the tensor 
(2.7) reads 

n 

pi') + N t x l (>> + p| s) ) = 0. (2.9) 

i=l 

Using the equations of state (2.8), we obtain from (2.9) the following integrals of motion 



= exp 



27o -J2 N * h< i 



const. (2.10) 



The Einstein equations R^ — RSfJ /2 = n 2 T^f (n 2 is the gravitational constant), can be written as 
R™ = k 2 [T^ - T6ff/(D - 2)]. Furthermore, we employ the equations flg - R/2 = k 2 T ° and R™> = 
k 2 [T™> - T6™;/(D - 2)]. Using (2.3)-(2.8), we obtain for them 

- Gi^iP + V = 0, (2.11) 
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'(70-7) = -K 2 Y, A(S 



,« E n k=1 N k hi SV 



s=l 



D-2 



x exp 



i=l 



Here 



Gij = N % 5 %3 - NiNj 



are the components of the minisuperspace metric, 



v = k 2 A(s) eX P 



]TiVA (s V-2( 7 - 7 o) 



(2.12) 
(2.13) 

(2.14) 



(2.10) is used to replace the densities in (2.11), (2.12) by expressions of the functions x l (t). 

After the gauge fixing 7 = Fix 1 , . . . , x n ) the equations of motion (2.12) are the Lagrange-Euler equations 
obtained from the Lagrangian 



(2.15) 



under the zero-energy constraint (2.12). 

Now we introduce an n-dimensional real vector space R n . By e%, . . . e n we denote the canonical basis in 
R° (ei — (1, 0, . . . , 0) etc.). Hereafter we use the following vectors: 
1. the vector x with components being the solution of the equations of motion 



x = x 1 {t)ei + . . . + x n (t)e r 
2. m vectors u s , each of them one component of the perfect fluid 



i=l 



D-2 



(«)' 



Let < ., . > be a symmetrical bilinear form defined on R n such that 

The form is nongenerated and the inverse matrix to (Gij) has the components 

6 ij 1 



Ni 2-D 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



The form < ., . > endows the space R n with a metric, the signature of which is (— , +, +) [p0[,pifl. Gij is 
used to introduce the covariant components of vectors u s 



(2.20) 



1=1 



Form them the bilinear form reads 



< u s ,u r >=Y^h!f ) hy ) N, 



i=l 



2-D 



(2.21) 



A vector y £ R n is called time-like, space-like or isotropic, if < y, y > is smaller, greater than or equal 
to zero, correspondingly. The vectors y and z are called orthogonal if < y,z >= 0. 
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Using the notation < ., . > and the vectors (2.16)-(2.17), we may write the zero-energy constraint (2.11) 
and the Lagrangian (2.15) in the form 



E 

2 



- < x,x > + K 2 e 2( ^ 7o) J2 A^e <u °' x> = 0, (2.22) 
2 s =i 

1 m 

L = -e 70 " 7 <x,x > -n 2 e~<-~< J2 A{s)e<U3,X> - ( 2 - 23 ) 



Furthermore, we take the so called harmonic time gauge, which implies 



7(*)=7b=X) JV * a:< - ^ 2 - 24 ) 

i=l 

^From the mathematical point of view the problem consist in solving the dynamical system, described 
by a Lagrangian of the general form 

1 m 

L=-<x,x>-Y J a {s) e <u *> x> , (2.25) 

2 8=1 

where x, u s G R n . It should be noted that the kinetic term < x,x > is not a positively definite bilinear 
form as it is usually the case in classical mechanics. Due to the pseudo-Euclidean signature (— , +, +) of 
the form < . , . > such systems may be called pseudo-Euclidean Toda-like systems as the potential given in 
(2.25) defines a well known in classical mechanics Toda lattices |35| . 

Note that, we have to integrate the equations of motion following from the Lagrangian (2.25) under the 
zero-energy constraint. Although an additional constant term — (with uq = £ R n ) in the Lagrangian 
(2.25) does not change the equations of motion, it nevertheless shifts the energy constraint from zero to 

1 m 

E= - < ±,x > +^a( s )exp[< u s ,x >] = -a( 0) = ~k 2 A^\ (2.26) 

s=l 

In our cosmological model, with (2^8) and ( 2.17 ), such a term corresponds to a perfect fluid with hf^ = 
for all i = 1, . . . , n. This is in fact just a Zeldovich (stiff) matter component, which can also be interpreted 
as a minimally coupled real scalar field. Taking into account the possible presence of Zeldovich matter, we 
have now to integrate the equations of motion for an arbitrary energy level E. 



3 Solving the equations of motion for the model 
which reduces to a classical open Toda chain 

We start from the Lagrangian (2.25) and the energy constraint (2.26) with 

n>m+l , m>2 . (3.1) 

The vectors u s may obey the relations 

<u s ,u s > = u 2 > , s = l,...,m , (3-2) 

<u r ,u r+ i> = -t^u 2 , ,r = l,...,m- 1 , (3.3) 

all the remaining <u r ,u s > = , (3-4) 

where u is an arbitrary non-zero real number. In this case the vectors u s are space-like, linearly independent, 
and can be interpreted as root vectors of the Lie algebra 
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A m = sl(m + 1, C). The Cartan matrix (K rs ) (see e.g. O, H) then reads 



(K rs ) = 



2 < u r , u s > 
< u r , u r > 



2 


-1 


. . 


. 





-1 


2 


-1 . . 


. 








-1 


2 . . 


. 











. . 


. 2 


-1 








. . 


. -1 


2 



(3.5) 



Now, we choose in R n a basis {/i, • • • , /«} with the following properties 

2m, 



s+l 



, s = 1, . . . ,m , 



with 



< u s , u s > 
< fi,fi> = Vn i i=l,...,n , 
</ s +i,/fe> = , < f k , fi >= rfki , s = l, ...,m ; k,l = m + 2, ...,n 

(rjij) =diag(-l, +!,...,+!) , i,j = l,...,n . 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



Note that, the basis contains, besides vectors proportional to u s , additional vectors / m +2 v ■ ■ iff n > m+1. 
By the decomposition 



/.r.t. this basis, with relations (3.2) - (3.4), ( |3.6| ) - ( |3.8| ) the Lagrangian (2.25) takes the form 

"m+l m 



1 / ,.W2 4 



-a«e 2 ^ 



£(<f) 2 -]>><f +1 

s=2 p=2 



fc=m+2 



r=3 



,(r-l) 2q r -q r - 1 -q r+1 _ (m) „2q m+1 -q" 



The equations of motion for g 1 (t), (7 m+2 (t),. . . ,q n {t) read 

g X (t) = 0, g m+2 = 0, q n (t) = 0- 

Then, 

q k (t) = a k t + b k , k = m + 2, . . . ,n. 
The other equations of motion for q 2 (t),. . . ,q m+1 (t) follow from the Lagrangian 



m+l 



s=2 p=2 

The linear transformation 



r=3 



>q s -lnC s , s = l,...,m, 

where the constants C ,. . . ,C m have to satisfy 

2^ X rs In C s = In — - — , r = l,...,m, 



(3.10) 



(3.11) 



(3.12) 



(3.13) 
(3.14) 



(3.15) 
(3.16) 
(3.17) 



brings the Lagrangian into the form 

m m — 1 



m — 1 



p=2 



(3.18) 
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The latter represents the Lagrangian of a Toda chain associated with the Lie algebra A m j35| when the 
root vectors are put into the Chevalley basis and coordinates describing the motion of the mass center are 

separated out. 

We use the method suggested in |jj for solving the equations of motion following from (|3. 18 ) and obtain 



m+l 



F s (t)= v ri ---v rs A 2 (ru---,r s )e^+-+^)t 



n<...<r 8 

where A 2 (ri, . . . ,r s ) denotes the square of the Vandermonde determinant 



A 2 (n,...,r s ) = ] [ (w rz 



The constants v r and w r have to satisfy the relations 

m+l 



H v r = A- 2 (l,...,m+l) , 



m+l 



w r = o 



The energy of the Toda chain described by this solution is given by 

m+l 

2 



m+l 



Finally, we obtain the folowing decomposition of the vector x(i) 

-2(lnF s (t)+lnC s ) 



c(t) = (a 1 t + b 1 )f 1 +J2' 



s=l 



< U s , U s > 



U S + Yl (aH + b k )f k 
k=m+2 



(3.19) 
(3.20) 

(3.21) 
(3.22) 

(3.23) 
(3.24) 



We remind the reader that the coordinates x l (t) of the vector x(t) are, with respect to the canonical basis 
in R n , the logarithms of the scale factors in the corresponding cosmological model. 
Let us introduce the vectors 



a = a 



7i+ E a fe /fc = E ale * eR " 

k=m+2 



i=l 



(3.25) 



(3.26) 



k=m+1 



with a 1 , (3 l being their coordinates with respect to the canonical basis. Using (3.7) and (3.8), we conclude 
these coordinates have to satisfy the following equations 



< a, u s >= ^ G lj a l u j ^ = , s = 1, . . . , m , 



< (3, u s >= Y G ijP iu \ s ) = > s = 1, . . . , m , 

where the u 1 ,^ are the coordinates of u s in the canonical basis (see (2.17)). 
The total energy E of the system is given by 



(3.27) 
(3.28) 



m+l 



E 



(3.29) 



G 



If m = n + 1, then < a, a >= - (a 1 ) 2 < 0. With fl3.29| ), we then obtain E < -^E . 

Finally, the scale factors of the multidimensional cosmological model with the Lagrangian (2.25) and the 
energy constraint (2.26) are given by 



m 

n [ p "(t) 



where 



F s (t) = C s ■ F s (t) , s = l,...,m. 
Using (2.10) we obtain the following solution for the densities of the perfect fluid components 

oW = A {1) e~ 2la — o {m) = A {m) c~ 2l ° Fm ~ 1 

P p2 ' P fi2 



P 



M =A M e -27o ViVi p=2,...,m-l. 



(3.30) 
(3.31) 



(3.32) 



where 70 is defined by (2.5) and may be calculated by (3.30). 

The constants C s are specified by (3.17). The solution contains the parameters a 1 , (3 l , v r , w r (i = 
l,...,n, r = l,...,m + 1) obeying the constraints (3. 27), (3.28), (3. 21), (3.22), (3. 29). If the energy E is 
arbitrary (see (2.26)) the solution has 2n free parameters as required. 



4 Example 

We consider a space-time manifold M with the following structure 

M = R x Ml x Mf x M 3 4 



(4.1) 



where the dimensions of Mf, Mf , and M3 are 3, 3, and 4, respectively. The first component of the perfect 
fluid shall have the Yip values 

(4.2) 
(4.3) 
(4.4) 



while the second component is given by 



,(2) 



h , h 



(2) 



4 X) = 



, 4 2) = o 



h is a real valued parameter with the restriction 



It is easy to check that the relations (^^), ( |3.3|) with s = 2 and r = 1 are fulfilled. 
In this case, the exact solution of the field equations gives the metric 



ds 2 = 



Fi(t)F 2 (t) 



xl - exp(8a i + 8(3 )dt 2 



ds\ ds\ 
P}(t) Pf(t) 



exp(2a t + 2f3 )dsl 



with the following definitions 



and 







3 




3 '^(t) 


F 2 (t) = k 2 




1 

3 




1 /i 2 F 2 (i) 



(4.5) 

(4.6) 
(4.7) 
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and 



F 1 {t) = Vl e Wlt +v 2 e W2t +v 3 e 

F 2 (t) = v 1 v 2 {w 1 - w 2 ) 2 e { ^ +w ^ + vtv^ - w 3 ) 2 e^ +w ^ 
+v 2 v 3 {w 2 -w 3 ) 2 e^ +w ^ . 



In our case, the energy E is given by 

E 



6a ° + zh? ^ + l " 2 + m 



(4.8) 
(4.9) 

(4.10) 



> means that Zeldovich matter is present in all the internal spaces (See the remarks at the end of 
sect. 3). With = ( 4.10D is the energy constraint specialized to our example. 

The nine parameters wi, w 2 , W3, v%, v 2 , U3, ao, /?o, E have to satisfy (4.9) and the two further relations 



and 



wi + w 2 + W3 = 



VlV 2 V 3 = [(wi - w 2 ){w 2 - w 3 ){wi - w 3 )}~ 



(See (|3.22|) and (|3.2l|)!). 

Finally, we have to give the expressions for the matter densities and p( 2 '. They read 



p (2) = A (2) 



-2/(3h)-2^pl-2/(3h)^ 



-2/(3h)-2, 



{t)F^ 2/{3h \t) 

p (2) _ A {2) 



and their quotient is 



,-8a t-8/3 



-8ao*-8/3o 



(1) 



AW F 2 {t) 



(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 



Although the solution is invariant with respect to the exchange (w±, w 2 , v\, v 2 ) — » {w 2 , w\, v 2 , v±) there is 
still enough freedom to build a lot of solutions that it is difficult to recognize many general properties of the 
solutions. What one can say is the following: We know that 



and 



V 3/l 


Fi{t) 


) K 


Fi(t) 



F 2 (t) 



F 2 {t) 



(4.16) 
(4.17) 



An easy but tedio us di scussi on of the different possibilities of choosing the parameters w\ and w 2 shows 
that the expressions ( 4.16 ) and ( 4.17 ) have for t — > ±00 the following shape: 



F 2 (t) 



and 



F 2 (t) 
Fi(t) 



Fi(t) 



t — *±oo 



t — *±oo 



f± oo (wi,w 2 )t 



g± a(w 1 ,w 2 )t 



(4.18) 



(4.19) 



with some functions /±oo, and g±<x> of the parameters wi and w 2 being negative for +00 and positive for 
— 00. This shows that the scale factors of the manifold Mi and M 2 go to zero for t — > ±00. Because of the 
behaviour of the scale factor of M\, as one could have expected, the solution can not properly describe the 
total cosmic evolution. But what one can do is arranging especially the parameters Wi and w 2 such that the 
behaviour of the scale factors is acceptable from a physical point of view (see below). 
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For the parameter values v\ = 1,V2 = 10, w\ = 1.1, W2 = 5 and two different values of h (hi = 2/3, hi = 2) 
(Zcldovich matter is present) the scale factors x l (t) are pictured on a certain time intervall. x\ is x l for 
h = hk- The integration constants are put equal to 1. 

As for the scale factor e x3 ^ of the manifold M 3 , we have 

(e* 3 ^) 3 oc |F 1 F 2 |e 3a °' lt (4.20) 

(the parameter f3 shall be absorbed by coordinates). For t — > ±oo we have 

|FiF 2 | -» ei*^.^)* (4.21) 

where k(w\,W2) is some positive function of the parameters wi and W2- 

The proper time T as a function of harmonic time t is given by integration of dT = e la dt with 70 = 
3a; 1 + 3a; 2 + 4a; 3 . For t — > ±00 the behaviour of both, the proper time T and the scale factor e x depends 
on the choice of a . This holds for the case that in all manifolds Zeldovich matter is present. If the case 
is excluded then ao is given by the energy constraint. We shall not discuss this point here because of the 
above mentioned restricted range of possibly physical importance of the solution. Instead, we present some 
drawings conveying how the scale factors and matter densities can behave for special values of wi, W2, v\, 
V2, and ao. 
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The energy densities pi and pi belonging to the parameter values v\ = 1, vi = 10, w\ = 1.1, W2 = 5 and 
two different values of h (hi = 2/3, hi = 2) (Zeldovich matter is present). p\ means the energy density pi 
belonging to hk- 
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